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CLOSE BINARY INTERACTIONS OF INTERMEDIATE-MASS BLACK HOLES: POSSIBLE ULTRA-LUMINOUS 
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ABSTRACT 

While many observed ultra-luminous X-ray sources (ULXs, Lx > lO^*^ erg s"') could be extragalactic X-ray 
binaries (XRBs) emitting close to the Eddington limit, the highest-luminosity ULXs (Lx > 3 x lO-'^ erg s~') 
exceed the isotropic Eddington luminosity for even high-stellar-mass accreting black hole XRBs. It has been 
suggested that these highest-luminosity ULXs may contain accreting intermediate-mass black hole (IMBH) 
binaries. We consider this hypothesis for dense, young 10^ yr) stellar clusters where we assume that a 50- 
500 M0 central IMBH has formed through runaway growth of a massive star. Using numerical simulations 
of the dynamics and evolution of the central black hole's captured companions, we obtain estimates of the 
incidence of mass transfer phases and possible ULX activity throughout the IMBH's evolutionary history. We 
find that, although it is common for the central black hole to acquire binary companions, there is a very low 
probability that these interacting binaries will become observable ULX sources. 
Subject headings: Stars: binaries: close — X-rays: binaries — Black Holes — Star Clusters 



L INTRODUCTION 

Extragalactic observations with the Einstein space X-ray 
telescope revealed a new category of sources, currently 
known as ultra-luminous X-ray sources (ULXs, Fabbiano 
1989). In 1990, surveys conducted with the ROSAT satel- 
lite were able to resolve many of these sources as dis- 
tinct objects not associated with emission from galactic nu- 
clei. Up to half of the galaxies surveyed were found 
to contain off-nuclear ULXs (ICo lbert & Mushotzkv 1999; 
[Roberts & Warwick' 2002; Lira, Lawrence, & Johnson 2000; 
IColbert & Ptak 2002; Colbert & Miller 2004). In analogy 
with X-ray sources observed in our own galaxy, many ULXs 
could be extragalactic X-ray binaries (XRBs) with stellar- 
mass black hole (BH) accretors (e.g., Rappaport, Podsiad- 
lowski, & Pfahl 2005). However, in many cases their lu- 
minosities exceed the Eddington limit for isotropic emission 
from an accreting stellar-mass BH. The lowest luminosity that 
commonly defines a ULX, Lx > 10^^ erg s~', corresponds to 
an Eddington black hole mass > 8M0 - still within the BH 
mass range allowed by stellar evolution models. Based on 
the X-ray survey conducted by Colbert & Ptak (2002), ULXs 
with Lx > 10-^^ erg s ~' are estimated to exist in ^ 12% of all 
galaxies ( Ptak & Colbert 2004) . Of the 87 X-ray (2-10 keV) 
sources observed, 45 had Lx > 3 x 10^^ erg s"'. This corre- 
sponds to the Eddington luminosity for the highest-mass black 
holes that can form through s ingle-star collapse, ^ 25Mpf 
for a metallicity ofZ = 0.001 (Belczvnski, Sadowski, & Rasio 
I2OO4.) . We see that as much as half of this observed ULX pop- 
ulation requires an alternative to the model of isotropic emis- 
sion from a stellar-mass XRB. 

Observations increasingly support the idea that some ULX 
sources may harbor intermediate-mass black holes (IMBHs, 
Colbert & Mushotzky 1999). Until recendy, no evi- 
dence existed for BHs in the wide mass range between 
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Stellar-mass and supermassive BHs. Many spectral and 
timing observations conducted with Chandra, ASCA, and 
XMM-Newton have confirmed that ULXs are not associ- 
ated with super novae and are c onsistent with accreting bi- 
nai-ies (.Fabbiano & White! l2003h . Additionally, ULX posi- 
tions are strongly correlated with young, star-forming re- 
gions. Examples include the nine ULXs found in the An- 
tennae ( Fabbiano, Zezas, & Murrav 2001) and the highest- 
luminosity ULX yet observed, in M82 (peak Lx 9 x 10"^" 
erg s"' , corresponding to an IMBH mass ~ 7OOM0 if Lx — 
L™; Kaaret et al. 2001; Matsumoto et al. 2001). The young, 
dense clusters where ULXs are often observed are also likely 
sites for IMBH formation ( Freitag. Giirkan. & Rasio 200jl 
Freita g. Rasio. & Baumgardi2005l: iGiirkan. Freitag. & Rasia 
■2004^ .' 

Accretion disk spectra provide a further test for the pres- 
ence of IMBHs: higher-mass BHs are expected to have cooler 
accre tion disks modeled as multi-color disks ( Mitsuda et a^ 
Some evidence for low-temperature disks was recently 
found using high-resolution XMM spectra of numerous ULXs, 
including those in NGC 1313, NGC 4559 Ho IX, and the 
Antennae ( Miller et alJl2003l l2004albh . As IColbert & Milled 
(2004) note, the sources observed to have cool disks are also 
the ones with luminosities too high to be explained by XRBs, 
even if they are beamed XRBs. Nevertheless, we note that 
X-ray spectral interpretation is subject to a number of model 
assumptions and may not be unique. 

Anisotropic, or beamed, disk emission is another possible 
explanation for apparently super-Eddington X-ray luminosi- 
ties. The degree to which emission is beamed can be defined 
by the beaming fraction b = VI /An, where Vt is the solid an- 
gle of emission. ULXs containing beamed stellar-mass XRBs 
could be a short phase of rapid mass transfer in the lifetime 
of ordinary XRBs, such as thermal-timescale mass transfer 
fKing et al. 2001). As demonstrated by King et al. (2001), by 
introducing mild beaming, the XRB falls within the stellar- 
mass BH range, assuming Lx < 10^" erg s~'. However, 
about 5-10% of the to tal ULX population is observed to have 
Lx > 10""^ erg s"' (Colbert & Miller '2004*: Pta k & ColberJ 
^04; Colbert & Ptak 2002). At these luminosities, very mas- 
sive BHs (~ 25Mq) could fall below the Eddington limit with 
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b < 0.25. To have the required mass ratio q = M2/M1 > 1, 
where Mi is the BH mass, these BHs would require a very 
massive companion. More moderate-mass BHs would require 
more severe beaming fractions that are not easily justified. 
Furthermore, in at least two instances observations seem to 
indicate isotropic X-ray emission. These include an isotropic 
nebula observed arond the ULX in Ho II (Pakull & Mironi 



120011) and the discovery of quasi-periodic oscillations (QPOs) 
from the brightest ULX in M82 ( Strohmaver & Mushotzkv 
l20()l . 

Transient behavior of a ULX source is an irnportant test for 
the presence of an IMBH ( Kalogera et al. 2004). If the ULX is 
a beamed stellar-mass XRB in thermal-timescale mass trans- 
fer, it is expected to produce persistent emission at or above 
the Eddington luminosity ( IKing et al.l f2001). A persistent 
IMBH ULX similarly would require a sustained M compa- 
rable to the Eddington mass transfer rate: 
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For the IMBH masses we consider, this is in excess of 
1O~^M0 yr~', corresponding to Lx ~ few x 10^" erg s~' . Tran- 
sient IMBH ULX sources are more plausible, because they 
confine the most extreme rates of mass transfer to short out- 
burst periods. Admittedly our current physical understand- 
ing of transient behavior is not complete. Nevertheless, it is 
generally accepted that an accreting binary becomes a tran- 
sient when the mass transfer rate driven by the donor is lower 
than a certain critical value and a thermal disk instability de- 
velops. In what follows we adopt the critical mass transfer 
rat e for transient beha vior of an X-ray uradiated disk given 
bv lDubus etaPdigggl) : 



M" 



1.5 X 10' 



Mi 

m;; 



-0.4 



R 



C 



10i"cm 
-0.5 



2.1 



5 X 10-4 



gs 



(2) 



for donor mass M2 and disk radius R. C is assumed to be 
5 X 10-4 for most donors. It can also be shown that the crit- 
ical M for transience corresponds to a minimum black hole 
mass that d epends o n companion mass an d orbital period 
jKing et al.il996.) . As lKalogera et al J (12004 have shown, this 
minimum BH mass is high for main sequence (MS) donor 
stars (1O^M0 on average for a lOM© donor). Red giant (RG) 
donors can more easily form transient systems with IMBHs. 
One should note that, due to the minimum mass limit, a 
stellar-mass BH binary is unlikely to be transient, so any tran- 
sient ULX observed is a strong IMBH candidate. 

In our simulations, we consider IMBHs that have 
formed through "runaway collisions" of MS stars. This 
is a channel for IMBH formation in which a series 
of rapid mergers occur as the cluster begins core col- 
lapse, causing growth of a central massive object that 
can collapse to a black hole (Port e^ies Zwart ^t ^ 1.1 
19991 fPortegies Zwart & McMilla nj l2002t 



Giirkan, Freitag, & Rasio 2004; Freitag. Rasio. & Baumgardtl 
2005: Freitag. Giirkan. & Rasio 2005). It has been sug- 
gested that an IMBH formed though runaway growth 
actually prevents core collapse and causes the core to 
reexpand ( Baumgardt, Makino, & Ebisuzaki 2004). In 
the dense clusters we are considering, the timescale 
for IMBH formation t hrough runaway mergers is esti- 
mated to be < 3 Myr (iPortegies Zwart & McMillanI l2002t 



iGiirkan. Freitag. & Rasiol 120041) . Because the core collapse 
time, tec, is expected to obey the relation tec ~ 0.15tie(0), 
where t,c(0) is the initial core relaxation time, this constrains 
tre to be < 30 Myr (Gurkan. Freitag. & Rasio 2004). 

Given that a plausible IMBH formation method exists, the 
question most pertinent to understanding ULXs is whether an 
IMBH, once formed, will gain close stellar companions that 
can sustain mass transfer at ULX X-ray luminosities. No prior 
studies have attempted to investigate numerically with both 
dynamics and binary evolution whether binary formation and 
mass transfer can occur with IMBHs in a dense cluster core 
(for an analytical exploration of tidal capture and its conse- 
quences, see Hopman, Portegies Zwart, & Alexander 2004). 
Here, we take the first step in studying IMBH binary inter- 
actions as a ULX possibility, using detailed numerical simu- 
lations that combine dynamical interactions with full binary 
stellar evolution. In § 2 we outline the details of our simu- 
lations. The results for our standard cluster model (details 
below) are presented in § 3. The variations of these results for 
other cluster models are described in §4. In § 5, we discuss 
our results, estimating a lower limit on the incidence of IMBH 
mass transfer, and we outline our goals for future simulations. 

2. CLUSTER SIMULATIONS WITH CENTRAL IMBH 

The goal of our simulations is to examine the ability of 
IMBHs in young, dense stellar clusters to form and main- 
tain close binary systems with companion stars. We em- 
ploy a numerical code developed by Ivanova et al. (2005, 
where a detailed description of the method can be found). 
This adopted hybrid code incorporates: (i) a stellar popula- 
tion synthesis code to evolve single stars and binary systems 
(StarTrack; Belczynski, Kalogera, & Bulik 2002; Belczynski 
et al. 2005 to be submitted); (ii) a semi-analytical prescrip- 
tion for 2-, 3-, and 4-body stellar collisions, disruptions, ex- 
changes, and tidal captures based on previously derived cross 
sections; and (iii) a numerical toolkit for direct N-body in- 
tegration of 3- and 4-body systems (Fewbody; Fregeau et al. 
2004). This combination provides us with a unique numeri- 
cal tool that lends itself to the study of binary interactions in 
cluster cores. At the young cluster ages we consider, large- 
scale cluster evolution has little effect, so certain simplifica- 
tions can be made. We employ a "fixed background" model, 
which holds the core stellar number density ne and velocity 
dispersion a constant (e.g.. Hut, McMillan, & Romani 1992 
Di Stefano & Rappaport 1994; Sigurdsson & Phinney 1995 
Portegies Zwart et al. 1997a; Portegies Zwart et al. 1997b 
Rasio et al. 2000). The cluster density profile is assigned 
by dividing the cluster into two regions: a dense core and an 
outer halo. This model is qualitatively based on mass segrega- 
tion. Note that the core is expected to expand on a very short 
timescale immediately after the IMBH is formed, smoothing 
out the cusp and supporting the assumption of a constant den- 
sity core (Baumgardt. Makin o. & Ebisuzaki 2004.) . 

In our models, an IMBH is added to the core at 1 Myr, in 
accordance with estimated IMBH formation times of 1-3 Myr 
in clusters (Portegies Zwart & McMillan 2002). To form an 
IMBH through runaway collisions, the initial core relaxation 
time, tie, must be < 25- 30 Myr (.Giirkan. Freitag. & Rasiol 
t). This time is given by: 



4.887rG2ln(7N)ne(m)2 ^ ^ 

dSnitzedl 19871) . where a^u is the three-dimensional velocity 
dispersion, ne is the core number density, N is the total num- 
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TABLE 1 
Cluster Models 



Model 


IMF 


tic [pc ^] 


A 


Kroupa 


1.33 X 10' 


B 


Flattened 


1.33 X 10' 


C 


Kroupa 


1.33 X 10"* 



Note. — The IMF in models B 
is a broken power-law mass func- 
tion 5(M) ~ M", with a = -1.25 
for 0.8 < M/Mq < 5 and o = -1.5 
for M/Mq > 5. Models A & B 
were tested with and without 3BBF. 
Model A is also tested with a range 
of IMBH masses from 50 to 500 
M0. 



ber of stars, (m) is the average stellar mass, and 7, the co- 
efficient in the Coulomb logarithm A = 7N, is assumed to 
be 0.01 (cf. Giirkan, Freitag, & Rasio 2004). For our stan- 
dard model (parameters are described below), trc — 29 Myr, 
in good agreement with the constraint for runaway collisions. 

We evolve the clusters to 100 Myr, as we consider only 
young clusters like those with which ULXs are most often 
associated. We can then neglect long-term contributions of 
the outer halo to cluster evolution and synthesize the entire 
stellar population of N = 2.7 x 10"* in the core. Mass segrega- 
tion within the core is incorporated through our choice of the 
initial mass function (IMF) used for the stellar population in 
the core. 

An initial binary fraction of 100% is assumed, which pro- 
vides the maximum possible concentration of binaries to in- 
teract with and become companions of the IMBH. Observa- 
tions indicate that young clusters are expected to have very 
high binary fraction s ( Levine, Lada, & Elston 1999; Apai 
|2P04; Delg ado-Dona te et al. 2004). This assumption is also 
consistent with the suggestion that cluster binary popula- 
tions a re depleted oyer time by interactions and stellar evolu- 
tion ( Mikkola'1983VHills'1984VBacon. Sigurdsso n. & DaviesI 
[1996; Fregeau et al. 2004; Ivanovaet al. 2005). Initial condi- 
tions for binary orbital parameters are assigned as outlined in 
§ 3.1 of .Ivanova et al.. ( .2005) . with slight modifications: 

• The binary mass ratio q = m2/m\ is assigned a uniform 
distribution < <7 < 1, where Mi and M2 are the binary 
primary and secondary, respectively. 

• A uniform logarithmic distribution is used for the bi- 
nary period, P, ranging from 0. 1 - 10^ d. 

• The thermal distribution is used for the binary eccen- 
tricities e with probability density p(e) = 2e. 

• Binary systems are rejected from the initial distribution 
if one of the binary components enters Roche lobe over- 
flow (RLOF) at pericenter. 



• We use the IMF ofKroupa' UMI) for primary stars be- 
tween 0.8 and 100 M0 in most simulations. 

• Secondary stars are restricted in a mass range between 
0.05 and lOOMg. 



Our standard cluster model (model A, see Table has a 1- 
D stellar velocity dispersion ct = 10 km s"' and a core stellar 
number density n^ = 1.33 x 10^ pc"^. The latter corresponds 
to a luminosity density > lO^L© pc"-' for a cluster age < 10*^ 
yr. The total cluster mass is ~ 5 x lO'^M©. We use 100 Mq 
as our standard IMBH mass, but we explore a range of IMBH 
masses from 50-500Mq (see §4.1). Note that each type of 
simulation must be repeated many times to obtain good statis- 
tics. 

We also consider two variations of our standard cluster 
model (see Table0: with models B we explore the effect of 
a flatter IMF, and with models C we examine results for less 
dense clusters. Detailed simulations of clusters that form an 
IMBH just before core collapse have shown that rapid mass 
segregation preceding IMBH formation creates a mass func- 
tion flatter than Kroupa's for stars > 5 - lOM© (Giirkan, Fre- 
itag, & Rasio 2004; Freitag, Giirkan, & Rasio 2005; Freitag 
2005, personal communication). These results have prompted 
us to consider the effects of mass segregation on the forma- 
tion of IMBH binaries. We have also examined whether any 
Brownian-like motion of the IMBH would cause it to inter- 
act with stellar populations that have different properties than 
those typical of the compact cluster cores. Assuming energy 
equipartition exists between the IMBH and the cluster, the 
IMBH in a constant-density core will experience radial oscil- 
lations of amplitude (Bahcall & Wolf 1976) 
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where /?core is the core radius, Mjmbh is the IMBH mass, and 
Mtot is the total mass of the rest of the cluster. For a 100 
Mq IMBH and total cluster mass ~ 5 x IO'^Mq, we obtain an 
oscillation amplitude that is safely small: R/Rmre — 2.5%. 

We do not normally allow binary formation through three- 
body interactions in any of our simulations, as the initial bi- 
nary fraction is 100%, but we do test the effect of 3-body bi- 
nary formation (3BBF) on models A and B. The details of our 
treatment of 3BBF are given in Iva nova et alJ ( 120051) . 

A large fraction of IMBH binaries acquire at least one, and 
often a series of tertiary objects; however, no population syn- 
thesis code currently exists that is capable of modeling full 
stellar evolution for triple systems. For triple systems that are 
considered stable after the detailed direct N-body modeling 
with Fewbody (see stability criterion in Mardling & Aarseth 
2001), we implement a "breakage" prescription that is phys- 
ically self-consistent. The energy required to eject the outer 
star is obtained by shrinking the inner binary orbit. The outer 
star is released unless the inner binary merges during shrink- 
age; in this case the inner system is allowed to merge and the 
outer companion is kept at its new, wider orbit to form the fi- 
nal binary system. The condition for merger is that the stars in 
the inner binary are in physical contact at pericenter or, in the 
case of an IMBH binary, the secondary star enters the black 
hole's tidal radius 



'Ml 

^tidal«/?2| ^ 



1/3 



(5) 



We recognize that this artificial treatment can affect the evo- 
lution of IMBH binaries, for example, by preventing eccen- 
tricity boosts via the Kozai mechanism ( Kozai 1962). For this 
reason we have undertaken a careful analysis of triple break- 
ing and its effect on our results (see § 3.5). It should also be 
noted that in some cases with a RG companion, release of the 
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0.4 0.6 
eccentricity 




Fig. 1. — Probability densities of IMBH binary properties throughout the 
simulations: companion mass, orbital separation and eccentricity. The ther- 
mal distiibution, p{e) = 2e, is shown on the eccentricity plot. Results from 
167 simulations of model A with a 100 Mq IMBH are shown. 



outer companion induces a common-envelope phase in the re- 
suhing binary. 

3. RESULTS FOR STANDARD MODEL 

In what follows we present the results of our standard, 
100 M0 IMBH cluster simulations (model A in Table 1) in 
the context of the formation of IMBH binaries that experi- 
ence mass transfer and may become observable as ULXs. In 
§3.1 we discuss the characteristics of IMBH companions, in- 
cluding their orbital parameters, evolutionary types, and mass 
transfer phases. In § 3.2 we present an example of a single 
simulation and more details on the history of the IMBH. In 
§ 3.3 we examine the X-ray luminosities produced by mass 
transfer and determine the time fraction over which the IMBH 
could be observable as a ULX. In § 3.4 we comment on the 
statistical interpretation of our results. Finally, in § 3.5 we an- 
alyze triple stellar systems and their possible effect on mass- 
transferring binaries. 

3.1. Statistical Results on IMBH Companions 

The distributions of three orbital parameters: companion 
mass, orbital separation, and eccentricity, are shown in Fig.n 
for all IOOM0 IMBH simulations for cluster model A. The 
prevalence of low-mass companions IMq) in the core is 
consistent with the model IMF that is dominated by such low- 
mass stars. Orbital separations of the dynamically-formed 
IMBH binaries range widely from ~ 1 - 10^ Rq. The lower 
limit corresponds to the Roche lobe radius of the IMBH com- 
panion. On the other end, the soft boundary for binaries is 
defined as the "escape radius" 
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In the simulations, a star is treated as a companion if its orbital 
separation is less than lOOasott, corresponding to the maxi- 
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Fig. 2. — Companion mass vs. orbital separation for 42 model A simula- 
tions with a 100 Mq IMBH, presented separately for the following compan- 
ion evolutionaiy types: MS, Hertzsprung gap (HG), first giant branch (RG), 
horizontal branch (HB), and asymptotic giant branch (AGB). Fewer He and 
WD companions undergo mass transfer, and they are not included here. Dots 
denote simulation timesteps when the IMBH has a binaiy companion. Open 
circles denote a mass-transferring companion (through RLOF). Open trian- 
gles denote binary companions associated with a triple system that was bro- 
ken. 



mum orbital separation, 10^ Rq. Fig.[2illustrates that the ma- 
jority of IMBH companions have orbital separations between 
100 and IQ*Rq. Eccentricities follow the initial thermal dis- 
tribution, which favors higher values (Fig.[0- These high ec- 
centricities tend to increase the occurrence of mass-transfer 
events in which the companion overflows its Roche lobe at 
pericenter 

We find that the IMBH has a companion for ~60% of the 
total simulation time, and it has a mass-transferring com- 
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panion for only ~3% of the total evolution time (Table |2}. 
Fig. 121 shows companion masses and orbital characteristics 
for different evolutionary types. The IMBH companions are 
overwhelmingly MS stars due to the much longer lifetime of 
this phase and their dominance by number. About 90% of 
unique 100 Mq companions are MS stars for at least part of 
their binary lifetime; the IMBH has an average of ^6 MS 
companions per simulation and only one post-MS compan- 
ion for every two simulations ["post-MS" will henceforth re- 
fer to all post-main-sequence evolutionary types capable of 
driving mass transfer - i.e., neutron stars (NSs) and BHs are 
excluded]. In contrast, post-MS companions are more likely 
candidates for RLOF than are MS companions. The percent- 
age of post-MS companions that undergo mass transfer in 
our simulations is much higher than the percentage of mass- 
transferring MS companions (see Tabled). As noted above, 
post-MS companions are also more likely to enter transient 
mass transfer phases, which are more conducive to ULX for- 
mation. However, because the MS lifetime is about 10 times 
longer, mass transfer phases of MS companions are typically 
of longer duration, and the numbers of MS and post-MS mass- 
transferring companions per simulation are comparable. (Ta- 
ble|3. 

3.2. Single Run Example 

The binary companion history of a single 100 M0 IMBH 
simulation, using model A, is shown in Figs.|3]&|4] This sim- 
ulation is atypical in that it has multiple mass transfer events, 
but the behavior of these events is characteristic of mass trans- 
fer seen in other simulations. Fig.l^shows the IMBH evolu- 
tion over time. The black hole has a companion for ~ 60% of 
its total evolution time, which is about the average we calcu- 
late for all 100 M0 simulations (Table|2j. Note the often high 
eccentricities of non-mass-transferring companions, which is 
also typical of our simulations (Fig.fQ. In the lower window 
of Fig. 13 one can see the hardening of several companion or- 
bits before they are disrupted by another interaction. The long 
mass transfer phase at the end of the simulation characterizes 
stable, MS mass transfer. The first companion that undergoes 
RLOF (at about 10^ yr) is the most interesting. Mass trans- 
fer begins on the MS and continues through the HG until the 
donor becomes a He HG star. After this, the star explodes 
and leaves a NS orbiting the IMBH. The companion's circu- 
larized orbit is kicked by the supernova to an eccentricity of 
about 0.6. As Fig. |3] shows, the NS is later forced into an 
even higher-e orbit by an encounter with an outer compan- 
ion shown in Fig.|4] This causes an eventual merger with the 
IMBH. The outer-companion swapping that occurs with the 
NS-IMBH binary is frequently-observed behavior for a close, 
stable binary; more details on the behavior of triple systems 
are presented in § 3.5. 

3.3. X-ray Luminosities 

We examine the range of X-ray luminosities produced by 
the various mass-transfer phases that occur in our standard 
model simulations. The binary evolution component of the 
simulations allows us to calculate the mass transfer rate M 
driven by RLOF donors to the IMBH. We compare this rate 
to the critical rate Merit for transient behavior (Eqn. |3 to de- 
termine whether the IMBH XRB is persistent or transient. If 
persistent, the transfer rate from the donor is converted to an 
expected X-ray luminosity; if transient, we assume that the 
XRB emits at its Eddington luminosity during disk outbursts 
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Fig. 3. — Evolutionary history of a single 100 Mq IMBH simulation. 
Top window: IMBH companion mass M2; middle window: eccentricity; bot- 
tom window: orbital separation. Circles denote timesteps with mass transfer; 
each non-mass-transferring companion appears at a different mass value. It 
is evident that the IMBH changes muliple binary companions during a single 
cluster simulation. Gaps in the top and bottom plots indicate times with no 
IMBH companion. 
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Fig. 4. — IMBH companion mass (M2) vs. orbital separation (a) for the 
same simulation shown in Fig. 1 . Dots denote timesteps with an IMBH com- 
panion; open circles denote mass-transferring companions. Triangles mark 
inner-binary orbital parameters of triple systems and solid lines connect to 
the outer orbital parameters. The maximum companion orbit radii for RLOF 
are shown at: the end of MS (line 1), the start of the RG branch (line 2), 
the maximum possible stellar radius (line 3), and the maximum stellar radius 
for clusters younger than 10^ yr (line 4). Line 5 marks the soft/hard binary 
boundary. 
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TABLE 2 
IMBH BINARY COMPANIONS 



MlMBH = 100 M0 



MC Runs 


167 


\^ 'L.omp / 


6.49 


(Nms> 


5.90 


\J-''postJVlS / 


0.47 


/mMT\ 
\"ms/ 


0.26 


/mMT \ 

V^postMS/ 


0.23 


% MSmt 


4.36 


% postMSMT 


50.00 


% tcomp 


57.63 


%tMS 


40.87 


% tpostMS 


4.44 


nz- fMT 


2.81 




0.11 



Note. — All simulations adopt cluster model 
A. MC Runs: number of distinct Monte Carlo sim- 
ulations evolved to 10** yr. (Ncomp): average num- 
ber of companions (of any type) per run. (Nms) & 
(NposiMs). average number of companions of each 
type per run. Most postMS companions are captured 
while on the MS and do not coirespond to differ- 
ent stellar companions. (Nj^|) & (N^^j^g): average 
number of MS and post-MS companions that undergo 
mass transfer per run. % MSmt & % postMSjyiT: 
number % of MS and post-MS companions that un- 
dergo mass transfer. % tcomp : average % of total simu- 
lation time spent with a companion of any type. % tms 
and % tposiMs : average % of time with MS and post- 
MS companions. % t^l & % t^^^^^^: average % of 
time spent in MS and post-MS mass transfer phases. 



(|Kalogera et al. 2004; Ivanova & Kalogera 2005). In model 
A clusters with a 100 Mq IMBH, the total mass-transfer time 
fraction is about 3%. Of this, about 10% is spent under tran- 
sient conditions, i.e., M < Merit- However, only a very small 
fraction of the transient mass transfer time, comparable to the 
transient duty cycle, will be spent in outburst. Observations 
of Galactic transients imply that this duty cycle is of the or- 
der of a few percent (McClintock & Remillard 2005). Such 
short duty cycles correspond to very low time fractions of the 
10** yr cluster age that could be associated with ULX emis- 
sion: < 0.03%. The probability of detecting one of these 
transient sources in outburst is unfortunately vanishingly low. 
The other 90% of mass transfer time in our simulations is per- 
sistent and a ULX will result, by definition, if the X-ray lu- 
minosity exceeds 10"*^ erg s"'. We find that a 100 M© IMBH 
becomes a persistent ULX for only ^2% of total mass trans- 
fer time in model A clusters; this corresponds to ~ 0.06% of 
the total age of the young clusters we consider here. We con- 
clude that either persistent or transient IMBH ULXs would 
be observable for only a tiny fraction (< 0.1%) of the cluster 
lifetime. 




fMT + tMT 
^MS ^ TpostMS 



Fig. 5. — Distribution of mass transfer time fraction for 167 standard model 
simulations with a 100 Mq IMBH. 



3.4. Averages vs. Medians 

It should be emphasized that the mass transfer time frac- 
tions we present are averaged over all simulations, and that 
individual simulations are given equal weight, but their re- 
sults do vary over a wide range. The majority of our simula- 
tions have no IMBH mass-transfer events at all; it turns out 
that the median mass-transfer time fraction for each cluster 
model is always at or very near zero (Fig.|5}. At the other ex- 
treme, IMBHs in a few model A simulations acquire very sta- 
ble mass-transferring companions and spend more than half 
of their total evolution time in mass transfer phases. It is clear 
that the averages reported are affected by the presence of these 
rare outliers, as one can see by comparing the entire model A 
data set (Table |2l tj^| and t^]^[jy[s) to the smaller set for which 
triple system data was obtained, which does not include the 
outliers (Table |3] tMx)- We would like to note that for the 
majority of our simulations the predictions for the IMBH bi- 
naries leading to mass transfer with X-ray luminosities in the 
ULX regime could be even lower than we report based on the 
average quantities and results. 

3.5. Triple Systems 

The significant number of triple systems that form and are 
artificially (although physically self-consistently, i.e., con- 
serving energy) disrupted in our simulations warrant a sepa- 
rate analysis of their possible effect on mass transfer The NS 
merger event shown in the example simulation (Figs.|3]and|3, 
which resulted from an eccentricity boost by a third object, in- 
dicates that three-body interactions can play an important role 
even when triples are artificially broken. By analyzing the ini- 
tial orbital characteristics of triples and their coincidence with 
mass-transfer phases, we develop a sense of whether our treat- 
ment of triples is, in fact, a reasonable approximation. 

Almost all of the triples are hard in the outer orbit (see 
Fig.0}, but depending on orbital parameters, the inner binaries 
may not be influenced by the outer companion. We quantify 
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TABLE 3 

Triple systems and IMBH companions 



M,, 



100 Mr. 



MC Runs 

(N,,,p/Nbi„> 
% tripciosc 
% MT,ff 
% MT,,o,e 
%tMT 



fMT 
'^aff 



* '■close 



109 
3.46 
0.50 

28.90 
47.72 
10.24 

1.60 
21.07 

4.96 



Note. — All simulations adopt cluster model A. 
MC Runs: number of distinct Monte Carlo simula- 
tions for which triples data were obtained (note that 
Tableslslandlslare based on smaller sample sizes than 
are other data). (Ntrfp): average number of unique 
triple systems per run. (N,rip/N(,in): average ratio 
of distinct triple systems to distinct binary systems 
formed. % tripciose- % of triples that are "close" 
(pout/uin < 3). % MTaff: % of MT companions 
that are potentially affected by one or more tertiary 
companions gained during their lifetime. % MT^osc : 
same fraction for close triples. % MT t^ff : percent of 
total MT time potentially affected by triples. % MT 
tciose^ same fraction for close triples. 



the degree of possible companion interaction using the ratio 
of the outer orbit's pericenter to the inner orbit's apocenter 
(pout/ain), that is, the closest possible approach between the 
two bodies. The result is shown in Fig.0(upper left window) 
and in Table |3] We assign pout/ain < 3 as the critical value 
of this ratio, below which it is likely that the outer companion 
will have a nonnegligible effect on the inner binary. These 
"close" systems comprise less than 1/3 of all triples. 

Another diagnostic for triple systems is to calculate the per- 
centage of mass transfer phases potentially affected. We con- 
sider "affected" mass transfer as that which follows or coin- 
cides with triple formation. More importantly, we also cal- 
culate the percentage of mass transfer phases associated with 
close triples. While almost half of mass-transferring compan- 
ions are affected, only ^ 10% are affected by close triples 
(Table|3}. If we instead consider the total mass transfer time 
that is similarly affected, this fraction drops to ^ 5%. 

Apart from calculating the fractions of binaries and mass 
transfer phases that could be affected by our treatment of 
triples, we also attempt to quantify the possible effects of 
triple breaking on the mass transfer phases. One way in which 
an outer companion may affect the inner binary is through the 
Kozai mechanism (K ozai 196 2) that can increase the eccen- 
tricity of the inner binary and lead to mass transfer. We can- 
not exclude a priori that by breaking the triples we miss some 
mass transfer events involving the IMBH. We calculate the 
timescale for the inner binary to reach maximum eccentricity 
due to the Kozai mechanism using the initial orbital parame- 
ters of the triples Jinnanen et alJIl 99*71 iMikkola & Tanikawal 



0.5 



0.5 



(1- 



(7) 



where / ~ 0.42 ln(l /eio)/(sin^(/o)-0.4)° ^ e^) and /q are the 
initial inner eccentricity and inclination, is the outer eccen- 
tricity. Mi is the total mass of the inner binary, Mq is the outer 
companion mass, ai and ao are the inner and outer orbital sep- 
arations, masses are in Mq, and distances are in AU. If t^oz is 
sufficiently short, the outer companion left alone could cause 
the inner binary to undergo mass transfer at pericenter How- 
ever, this would be possible only if r^oz is shorter than the 
colhsion timescale, t^oh, appropriate for the triple system: 



TcoU ■ 



:3.4xlO'^xP//Xfn^^ri 

I 



1+913 



M,o,+ (M) 



"■^d "^tot '^lOJ 



(8) 



where (M) is the average single star mass, is the clus- 
ter's stellar number density in units of lO^pc"-', vio = 
Voo/(10km s"'), v'oo is the relative velocity at infinity, Pd is 
the binary period in days, and a strong encou nter is defined 
by the closest approach d^ax < ka, with ^ ~ 2 flvanova et alJ 
I2005h . We find that only ~ 1 1% of all triples have Tkoz < Tcoii- 
If these systems were not artificially disrupted, they could po- 
tentially induce mass transfer through Kozai resonance that 
is not seen in our simulations. However, we also note that 
disrupting a triple changes the properties of the inner binary, 
and this can also lead to mass transfer. We find that about 
7% of all tertiary companions are broken before or during 
a mass tranfer phase in inner binaries; this is the maximum 
fraction of broken triples in our set of standard simulations 
that could artificially create mass transfer by shrinking the in- 
ner binary. Furthermore, we have examined the simulation 
data and found that there is essentially no overlap between 
disrupted triples that precede mass transfer and triples with 
Tkoz < T'coii. Because the two populations are comparable in 
size (7% and 11%, respectively), we conclude that these pos- 
sible effects on mass transfer events more or less mutually 
cancel. 

It should be noted that a small fraction 4- 8%) of triple 
systems form when a tight binary is captured by the IMBH; in 
this case, the IMBH is recorded as the outer, "ejected" com- 
panion and binary evolution is subsequently followed for only 
the two closely orbiting companions. The data from this small 
number of unusual disrupted triples are not included in our re- 
sults and has a negligible effect on the overall IMBH binary 
evolution. 

4. PARAMETER STUDY OF CLUSTERS WITH CENTRAL 
IMBH 

4.1. IMBH Mass 

Table 0] outlines the evolutionary differences between 
IMBH masses in model A clusters. The overall increase in 
companion capture with IMBH mass is apparent. Note that 
the fractions of both MS and post-MS companions that un- 
dergo mass transfer are fairly consistent for all IMBH masses. 
This underscores the greatly increased likelihood for a post- 
MS companion to undergo mass transfer compared to a MS 
companion, independent of external factors. The time frac- 
tions spent with a companion of any type also increase mono- 
tonically with IMBH mass, but because BH companions are 
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TABLE 4 
IMBH BINARY COMPANIONS 



MiMBH [M©] 


50 


100 


150 


200 


500 


MP Runs 




iO / 


7Q 




44 


/M \ 
\-L>comp / 




D.4y 


1 n ^7 
iu. JZ 




ZZ.OD 




2 07 








20 98 


U^postMS/ 


U. 






U.J J 






0.10 


0.26 


0.46 


0.48 


0.91 




0.14 


0.23 


0.38 


0.21 


0.32 


% MSmt 


5.00 


4.36 


4.76 


4.25 


4.33 


% postMSMT 


40.00 


50.00 


45.45 


38.89 


32.56 


% tcomp 


24.98 


57.63 


70.04 


76.51 


89.69 


%tMS 


18.76 


40.87 


53.69 


60.89 


56.95 


% tpostMS 


3.59 


4.44 


4.05 


1.77 


2.36 


% t'^T 

* 'ms 


0.21 


2.81 


0.76 


1.71 


2.43 


DOstMS 


0.02 


0.11 


0.16 


0.05 


0.10 



Note. — All simulations adopt cluster model A. Abbre- 
viations given in Tablel2l 
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Fig. 6. — Increase in companion acquiring rate for a typical 500 M© sim- 
ulation over a 100 M© simulation (cf. Fig. 1) is shown; the IMBH has 
a companion for almost all of its lifetime. Top window: companion mass; 
middle window: eccentricity; bottom window: orbital separation. 



more prevalent at higher masses, this behavior is not as clearly 
seen for MS or post-MS companions. The evolutionary his- 
tory of a typical single 500 Mq simulation is shown in Fig.|S] 
to illustrate the large fraction of its lifetime spent with a com- 
panion. 

A clear result in Table|4]is that 50 M© IMBHs have a very 
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Fig. 7. — Triple systems formed for 16 simulations of 100, 150, 200, and 
500 Mq IMBHs. The outer orbit's eccentricity is plotted vs. the ratio of 
orbital separations for each triple system. Systems for which Pom/fin is < 3 
are denoted by open triangles; other triples are closed circles. 



low rate of companion capture and mass transfer events. This 
places a lower limit on the IMBH mass likely to form a lu- 
minous X-ray binary. For higher IMBH masses, we empha- 
size as in § 3.4 that the average mass transfer time fractions 
are typically dominated by a few simulations with very long 
mass transfer phases. Consequently, we do not observe any 
significant difference in f^J or f^^JtMS between 100-500 M© 
IMBHs. 

As in 100 M© IMBH clusters. X-ray luminosities for other 
IMBH masses produce transient or persistent ULXs for in- 
significant fractions of the cluster lifetime, except possibly for 
the 500 M©. Clusters with 500 M© IMBHs do spend signifi- 
cantly more of their mass transfer time as ULX sources; these 
systems are transient for ~40% of mass transfer time and per- 
sistent ULXs for ~12% of mass transfer time. Even so, these 
ULXs would be observable for only ~ 0.35% of the total clus- 
ter lifetime. 

The rate of triple system formation increases monotonically 
with IMBH mass (Fig.0 Table g). Clusters with a 50 M© 
IMBH have very few triple systems, as is expected based on 
their low rate of IMBH binary formation. 500 M© IMBHs 
have by far the highest average rate of triple formation. In 
fact, more tertiary systems than binary systems form in this 
case, in a ratio of L6 : 1 as compared to 0.5 : 1 for 100 M© 
simulations (Table |5}- Our artificial treatment of triples limits 
our ability to analyze reliably clusters with an IMBH >500 
M©. In the 100-200 M© range, we do not find any conclusive 
pattern in the percent of triples with close orbits, or in the 
percent of mass transfer affected by triples. 

4.2. Cluster Models 

To determine the dependence of our results on initial clus- 
ter conditions, we test two variations of our standard cluster 
model, as outlined in Table ^ We use our standard IMBH 
mass of 100 M© for all these simulations. 
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TABLES 

Triple systems and IMBH companions 



TABLE 6 

Binary companions of 100 M0 IMBHs 



IMBH Mass [Mq] 


50 


100 


150 


200 


500 


MC Runs 


20 


109 


53 


16 


39 


(N,„p) 


0.90 


3.46 


7.98 


9.44 


35.59 


(Ninp/Nbin> 


0.38 


0.50 


0.75 


0.81 


1.58 


% tripciosc 


50.00 


28.90 


18.91 


26.49 


22.05 


% MTrff 


66.67 


47.72 


62.22 


33.33 


71.43 


% MT,to,e 


0.00 


10.24 


13.33 


6.67 


20.41 


%tMT 


0.32 


1.60 


1.11 


1.00 


2.02 


% t'^T 

* W 


3.13 


21.07 


41.82 


25.16 


67.08 


or tMT 


0.00 


4.96 


7.57 


0.00 


27.14 



Note. — All simulations adopt cluster model A. Abbreviations 
given in Tablel3l 




Fig. 8. — Probability density of companion mass at each timestep, shown 
for cluster models A (167 runs), B (102 runs), A+ (12 runs), and B+ (19 runs). 



In model C, when a lower number density is used (n^ =10'* 
pc"-'), the IMBH has a companion for only ~ 35% of its evo- 
lution time, about half the companion time fraction calculated 
for the higher-density model (Table |6j. The IMBH spends 
only a fraction of a percent of total simulation time with mass- 
transferring companions (0.53%). We conclude that clusters 
with our standard core number density, nc = lO*" pc"^, are 
much more likely to contain an accreting-IMBH ULX. 

IMBH companions in clusters with a flattened IMF (model 
B) are slightly more massive, as expected, with most between 
^ 1 - IOM0 (Fig. |8}. The IMBH has a binary companion 
for a shorter time fraction (40%) than in model A simulations 
(Table |6j. A decrease is also seen in the mass transfer time 
fraction between the two models. This could possibly result 



Model 


A 


A+ 


B 


B+ 


C 


MC Runs 


167 


12 


102 


19 


37 


(Ncomp) 


6.49 


7.75 


4.35 


4.32 


2.16 


(Nms) 


5.90 


6.92 


3.15 


3.05 


2.03 


(NpostMs) 


0.47 


0.50 


0.59 


0.58 


0.13 




0.26 


0.50 


0.16 


0.21 


0.05 


V^postMS/ 


0.23 


0.25 


0.21 


0.21 


0.03 


% MSmt 


4.36 


7.23 


4.97 


6.90 


2.67 


% postMSiviT 


50.00 


50.00 


35.00 


36.36 


20.00 


% tcomp 


57.63 


61.43 


40.35 


39.82 


34.71 




40.87 


48.90 


18.88 


23.19 


32.35 


% tpostMS 


4.44 


5.56 


4.69 


2.10 


0.50 


Of (MT 


2.81 


1.43 


0.54 


0.17 


0.52 


* 'postMS 


0.11 


0.05 


0.04 


0.07 


0.01 



Note. — See Table^for details of each model. Mod- 
els A+ and B+ are identical to A and B, but with 3BBF 
allowed. Abbreviations given in Tablel2l 



from having more massive stars in the core, due to the flatter 
IMF, that interact with and can disrupt IMBH binaries. How- 
ever, we caution again that the mass transfer time fraction for 
our standard model is dominated by a few outliers, so that the 
difference seen in model B may or may not be significant. 

The time in which the model B clusters are observable as 
ULX sources is very short, as in model A clusters: a transient 
ULX forms for about 3.5% of total mass transfer time, and 
persistent mass transfer at ULX luminosities occurs for about 
7% of mass transfer time. These imply a possible ULX time 
fraction of < 0.05%. Furthermore, the behavior of triple sys- 
tems in model B is qualitatively and quantitatively similar to 
that of triple systems in standard model clusters. 

Small differences can be seen when 3BBF is allowed in the 
cluster 3BBF enhances the rate at which massive objects can 
grow through collisions, since the hard formed binaries can 
collide with other formed binaries. As a result, companions 
are a bit more massive in clusters with 3BBF (Fig. |8}. The 
difference in the mass distributions is small, but a larger ef- 
fect is seen in the mass transfer time fraction. The time the 
IMBH spends with a mass-transferring companion decreases 
by more than half when 3BBF is included. This phenomenon 
is probably similar to that which causes the disparity in mass 
transfer time between model A & B clusters. 

Due to the 100% initial binary fraction, we do not need 
3BBF to form binaries. Overall, we see that our cluster mod- 
els are qualitatively similar with and without 3BBF, so we 
conclude that we can safely omit 3BBF from our simulations. 

5. DISCUSSION 

We have undertaken the first detailed investigation of IMBH 
interactions in young (~ 100 Myr), dense cluster cores, in- 
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corporating binary dynamics with full binary star evolution. 
Specifically, we study the prevalence and characteristics of 
interacting binaries involving the central IMBH, and we ex- 
amine and can better understand their role as potential ULX 
sources. We use numerical cluster simulations that simplify 
some aspects of detailed dynamical evolution and allow us 
to focus on individual core stellar interactions while largely 
ignoring long-term cluster evolution. We analyze our re- 
sults in terms of the frequency of IMBH binaries formed and 
their properties, the frequency of associated mass-transferring 
IMBH binaries, and the associated X-ray luminosities ex- 
pected. With our cluster model A, we explore an IMBH mass 
range of 50-500 M0, and we study two variations on this 
model as well (Table [Q. Finally, we analyze the nature of 
triple systems formed as a possible influence on our results 
here and as a guide for designing future simulations. 
Our main results are summarized as follows: 

• IMBHs more massive than 100 M0 in model A clusters 
spend the majority (> 50%) of their evolution time with 
a binary companion. 

• Mass-transferring IMBH companions, though, are rel- 
atively rare (~ 3% of the IMBH evolution time) com- 
pared to the total time the IMBH spends with a com- 
panion. 

• IMBHs spend more time in mass transfer phases with 
MS companions than with post-MS companions. 

• A much larger number fraction of post-MS companions 
than MS companions undergo mass transfer. 

• In model A clusters, IMBH binaries are found to be per- 
sistent X-ray sources for about 90% of the total mass 
transfer time, but the X-ray luminosity reaches ULX 
values (Lx > 10^'' ergs"') for only ~2% of the mass- 
transfer time. For the rest of the mass-transfer time 
(~10%), IMBH binaries are found to be transients, 
which have short duty cycles. Given that the typi- 
cal total mass-transfer time is ~ 3% of IMBH evolu- 
tion time, these results amount to a vanishingly small 
observable-ULX time fraction (< 1%) over the lifetime 
of the young clusters we consider 

• Higher-mass IMBHs capture more companions 
throughout their evolution, making them more likely 
in a statistical sense to capture a companion that can 
overflow its Roche lobe. 

• Model A clusters with a 5OM0 IMBH and cluster 
model C have very few mass-transferring IMBH bina- 
ries. 

• In model B clusters, the IMBH spends less time with 
a companion (40% of total time) and less time in mass 
transfer (< 1% of total time) than in Model A clusters. 

• For all cluster models and IMBH masses explored here, 
the possible ULX time fraction of the cluster lifetime 
(lOOMyr) is found to be typically of 0.1% or even 
smaller 

We note that results prese nted by 

iHopman. Portegies Zwart. & A lexander' (2004) imply a 
higher incidence of possible ULX activity associated with 
MS companions to central IMBHs (typically with duration 



of > 10^ yr). However, these calculations are based on 
analytical estimates of the possible tidal capture rates for an 
IMBH in a cluster of single MS stars and the subsequent 
isolated binary evolution once a MS companion is acquired. 
These binary evolution calculations are not integrated with 
any dynamical evolution (analytically or numerically) and, 
based on our results, we conclude that the assumption of 
isolated evolution is not realistic. Given the substantial 
differences in calculation methods, however, it is not possible 
to make direct and quantitative c omparisons between the 
study bv lHopman. Portegies Zwa rE & Alexander (2004) and 
our study. We note that in our simulations we do allow for 
tidal captures (for details of implementation see Ivanova 
2005), but they turn out to be a negligible contributor to the 
rate at which IMBHs acquire companions (exchanges into 
binaries is by far the dominant process). 

Last, we recognize triple systems as one possibly limiting 
factor in our simulations. Since triples are artificially dis- 
rupted as they form, we examine the influence of triples in 
terms of their initial orbital parameters and their coincidence 
with mass transfer phases. We find that while a large num- 
ber of triple systems form per simulation, and while about 
1/3 - 2/3 of all 100-200 Mq mass-transferring companions 
do coincide with triple system formation, only a small fraction 
of mass transfer phases are affected by triples with "close" or- 
bits. Furthermore, two opposing phenomena related to triple 
breaking may cause some uncertainties. Mass transfer may 
be prevented when potential Kozai evolution is artificially in- 
terrupted, or it may be artificially induced when release of the 
outer companion shrinks the inner binary. We find that both 
effects are relevant for only small, roughly equivalent frac- 
tions of triples, causing a negligible net impact on our mass 
transfer results. Therefore, we conclude that our treatment of 
triples in the regime of IMBH masses considered here does 
not influence our conclusions about the possibility of IMBH 
ULX formation. However, given the higher occurrence of 
triples in our 500 Mq runs, the issues of multiples and cusp 
formation are most probably more important for more mas- 
sive IMBHs (> 500 Mq) than those considered here. 

For these more massive IMBHs, a more physical treatement 
of multiple systems would in principle allow the remaining 
uncertainties regarding triples, multiples, and the ultimate for- 
mation of central cusps to be resolved. However, to examine 
ULX formation, treatment of the stellar evolution of multi- 
ple systems cannot be ignored (especially when it involves 
massive-star companions that lose mass in winds, experience 
core-collapse and can themselves disrupt orbits and alter dy- 
namical evolution). Simulations that incorporate all of these 
effects are simply not possible at present. Instead, for our 
future explorations of the IMBH-ULX connection in dense 
clusters, we plan to implement simulations that would repre- 
sent an improvement over the ones presented here. Such fu- 
ture calculations could differentiate between triple systems in 
which the inner binary may be disturbed and those in which 
the outer companion may be ignored. In the latter case, the 
third companion would be kept in the system and its dynami- 
cal evolution alone could be simulated, while full stellar evo- 
lution and dynamics are calculated for the inner binary. 
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